Time-evolution stability of order parameters and phase diagrams of bosons on optical 
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Stemming from the Heisenberg equations of motion, we study the time-evolution stability of the 
order parameters for the cold atoms on optical lattices. The requirement of this stability of the order 
parameters endows the phase diagram with a fruitful structure in the superfluid phase. For the one- 
component Bose-Hubbard model, we see that this stability of order parameter leads to a physically 
receivable phase diagram. For two-component bosons, we show that the molecules are preformed in 
the atomic superfluid and then condenses into a molecular superfluid phase at a critical repulsive 
inter-species interaction, which resembles the pre-pairing mechanism in high T c superconductor of 
Cu-O cuprates. 
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Ultra-cold atoms on optical lattices have offered a 
paradigmatic scene to study strongly correlated phe- 
nomenon in a highly controllable environment 00 y El- 
Superfluid/Mott-insulator transition was witnessed on 
87 Rb atoms in a three- and one-dimensional lattice 
0. Theoretically, the strongly correlated phenomenon 
for boson systems may be studied by the Bose-Hubbard 
model 0, and experimental feasibility was also estab- 
lished by microscopic calculations of the model parame- 
ters for cold boson atoms in optical lattices [|J. 

Very fruitful structures of the phase diagram have been 
explored . The phase boundary between the superfluid 
and the Mott insulator has been determined and ver- 
ified by the experiments 0, 0- In the Mott insulator, 
the Bose-Hubbard model can be mapped into an effec- 
tive spin Hamiltonian and the phase structure has been 
extensively studied Q. On the superfluid side, it is also 
predicted that there are complex quantum phase transi- 
tions hei 

The Bose-Hubbard model can only be exactly solved 
in one-dimension for a limited case. Various approxi- 
mations have been applied. If an approximation ground 
state describes a spontaneous breaking of symmetry, it 
has to be stable in the time-evolution. A powerful tool 
to study this stability is the Heisenberg equations of mo- 
tion of the corresponding order parameter operator. For 
any equilibrium state described by a time-independent 
density matrix p with [H, p] = 0, the Heisenberg equa- 
tion of motion of an operator O yields 



i(0) = ([H,0]) =Tr(p[H,O])=0, 



(1) 



which gives the exact constraint on order parameter: 
When an approximation state is readily to describe an 
ordered state, one has to check that the corresponding 
order parameter is stable or not in the sense obeying Q . 
The first such example has been provided by Zhang for 
the Hubbard model of the electrons and it has been gen- 
eralized to the t-J model These works suggested 
that the Hubbard model may not be a good model for 
the high T c superconductor while the t-J model is possi- 
ble. 



For bosons on the lattice, we shall show that similar 
exact constraints on the order parameters play an 
important role in the study of the phase structures of 
the systems. We shall use the variational wave function 
which yields the mean field state for Bose-Hubbard 
systems 0, ^| . It is seen that the variational energy 
will be affected by the constraint on the corresponding 
order parameter. We shall see that even for the simplest 
one-component Bose-Hubbard model, these exact results 
lead to a more physically receivable phase diagram. 
For the two-component Bose-Hubbard model with a 
commensurate filling, the most intriguing finding is that 
except the up-critical point from the superfluid to the 
Mott insulator, there is a phase transition from the 
atomic superfluid (AS) to the molecular superfluid(MS) 
in a low-critical point for a repulsive inter-species 
interaction. The existence of the MS phase for U > 
is because the local rij = 1 is violated in the superfluid 
and V > U favors two atoms in different species to form 
a molecule. Our analysis will show that the molecules 
are formed in the AS phase and condensed as U reduces 
to a critical value, which resembles the prepairing 
mechanism for the high T c superconductivity. To be 
specific, we confine our discussions to the single atom 
filling per lattice site in this work. However, it is easy 
to be generalized to other integer filling factors and the 
mixture of the Bose- Fermi atoms. 

One- component model: We now begin with the one- 
component Bose-Hubbard model, which is described by 
the Hamiltonian 

H x = -ty^(a\aj +h.c.) + X! n * - ^ - M^ n,; '( 2 ) 

(ij) i i 

where the notation is standard. For the uniform order 
parameters, the Heisenberg equations of motion for a and 
a 2 give 

— i(a) = tz(a) — V{na) + p(a), 

-i(a 2 ) = 2tz(a s a) + 2p(a 2 ) - V(na 2 }, (3) 
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where as is the boson operator at a nearest neigh- 
bor and z is the partition number. Precisely, the su- 
perfluid and dimerized order parameters are defined by 
(a) = Tr(pa) and (a 2 ) = Tr(/?a 2 ). In the zero tempera- 
ture limit, they are reduced to the expectation value cor- 
responding to the ground state. In the dilute gas limit, 
a widely accepted approxiamtion is, instead of the ex- 
act ground state, to introduce a site factorizable wave 
function which takes the form 03 



i*i>=n 



"it "if. Xi 

sin — aj + cos — sin — 
2 1 2 I 2 



1 Xi |2 
—= cos —a\ 
V2 2 1 



10) 



We confine the filling number per site to be one through- 
out this work. The order parameters now are (a) = 
(2 + V2) sin 0/(2^2) and (a 2 ) = cos 2 I /V2 since \ = n/2 
when the order parameters do not vanish and (n) = 1. If 
we do not consider the time evolution of the order param- 
eters, there are three possible phase: the Mott insulator 
phase with 9 = tt ((a) = (a 2 ) = 0), the dimerized phase 
with 9 = ((a) = and (a 2 ) 7^ 0) and the superfluid 
phase with both (a) 7^ and (a 2 ) 7^ 0. If taking the 
chemical potential fi as a constant, the previous varia- 
tional result has shown the critical points of the SF/MI 
and SF/D are in V/zt w ±5.83 7]. The plus one is 
reasonable and has been examined experimentally. The 
minus one seems to be a mean field artifact because it was 
known that in the one-dimensional model, the dimerized 
state dominates when V/zt < — 1 for V < To 
improve the result from this simple mean field theory, 
we consider the stability of the order parameters in time 
evolution. There are now two different variational states 
in superfluid phase. One is that (a) is stable and (a 2 ) 
is not and another is vice verse. The former subjects 
to the constraint —i(a) = which turns the first equa- 
tion of |0| to ^ = fii = -tz + V(y/2 - 1). The latter to 
—i(a 2 ) — 0, which turns to /1 = [12 = 
The variational energies to be minimized are 



where 0(x) is the step function. We depict these energies 
in Fig. 1(a). It is found that in the whole range where 
(a) 7^ 0, E 2 > Ei. However, at the point V/zt w —0.4, 
E 2 = Ei. This gives that the lower critical point of 
SF/D transition. The whole phase diagram of the system 
is shown in Fig. 1(b). There are three phases which 
are identical to those in Fig. 1 (b). The phase with 
stable (a 2 ) 7^ and nonzero but unstable (a) does not 
exist. However, due to the energy of the dimerized phase 
is lowered, the phase boundary between SF-D and D is 
(4) shifted a much reasonable value V/tz ss —0.4. 

In our variational state, the states with three or more 
atoms per site have been neglected. For a repulsive 
interaction, this approximation is appropriate because 
there was a large three-body collision loss in experiment. 
For a large negative V or denser atoms, the order 
parameters (a n ) (n > 3) are possibly stable. We can use 
the procedure provided before to discuss. Our procedure 
can also be applied to the other integer filling by taking 
an appreciate variational wave function. 

Two- component model: Turn to the two-component 
Bose-Hubbard model. The Hamiltonian is given by 

H 2 = - ^ tcr( a la a jcy + h.c.) + - ^ (n io - - 1) 

+ J7^(n <T - \){nn -h-^2^n icr . (8) 



E sl M0) = -^P^tz sinK - - 
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Mi,2- (5) 



Minimizing these variational energies, the optimizing pa- 
rameter is given by 



cosfli 



V 



tz{3 + 2V2)'' 



COS 02 = V2 



V 



tz(3 + 2V2) 



.(6) 



The requirement | cos ^1.2 1 < 1 gives the regimes of 
the minimal Ei S)2s defined: V^ ] < V < V$ with 

V^/zt w ±5.83 and V^ ±] /zt « 5.83( v / 2±l). Outside of 
these ranges, the energies are defined by Ei^(0 = 0,n). 
Thus, corresponding to two cases, the optimizing energies 
are 

£1,2 = s sl , s2 (0 = o)©(^2 ) -y) 

+ E shs2 {V/tz)®{V - V$)Q{V$ - V) 

+ E s i, s2 (6 = 7r)e(V-Vi { p), (7) 



Generally, the tunneling amplitudes t^i, the intra-species 
interactions V\ t \ may be different. It was known that if 
V\Vi < U 2 , one can not minimize the ground state energy 
with two nonzero occupation numbers (n\ \). The system 
must be phase separated pjj. We shall not discuss the 
phase separation here and consider only the interactions 
restricted in V a > and V^Vy > U 2 if U > 0. For U < 0, 
the phase diagram has been extensively studied and we 
shall not go to the details here 0, ^| . The total atom 
number is conserved in the system and the ratio between 
the different species atom numbers is self-consistently de- 
termined by the internal field and the constraints given 
rise to by the Heisenberg equations of motion. We now 
have five possible order parameters, (a CT ), (afaj.) and (a 2 ) 
if we fix the filling factor per site to one and neglect the 
three or more atom occupation in a single site. The time- 
evolutions of the order parameters are given by 

—i{a a ) = (zt a + fi a + —){a a ) - U {n s a a ) - Vrin^a^), 

~i(a 2 a ) = (2zt a + 2na + U){a s , a a a ) - 2C/(n ff a 2 ) 
- V a {n a a 2 a ), 

-«((aT a l)> = \ z (h + < l)(( a <5,T a l> + ( a <5.1 a T>) 
V 

+ ( 2" ~ U )( na 1 a l) + (MT + Ml)( a T a i)- ( 9 ) 
The most general factorable variational wave function in 
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our case is that 



And the variational energy reads 



l$2 
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, . Xi t , Xi t 

" m ^[ sin y a lt + cos y a ii 



1 <os ~2 l sm T s ~2 l sm "2 a »T a a 



1 / . Ci +2 Ci t2 

+ — ^ I sin — a'- t + cos —a 



V2\ 2 ^ 



|0). (10) 



Note that this factorable assumption implies that an ap- 
proximation neglecting all the spatial correlations. Thus, 
the BCS-type state may not appear in our phase diagram 
although it is a possible ground state Q . The condition 
(ri) = 1 fixes r\ = 7r/2 if ^ it. The Mott insulator phase 
for this Hamiltonian and the variational wave function 
have been extensively studied by authors of Ref. . We 
here focus on the homogenous state although we can not 
exclude the inhomogenous states such as spin glass and so 
no. They are the atomic superfluid state with (a a ) =/= 0, 
the molecular superfluid state with (a^a^) ^ and the 
dimerized state with (a 2 ) ^ 0. 

In terms of ©, the constraint identities may be ob- 
tained by requiring the stability of these order parame- 
ters. The variational energy now is given by 



-i T z(a T ) 2 - t lZ ( ai } 2 + U (n T n x ) - — 
(nf) + (nfj 



(11) 



(Vr/2 + MT )(n T ) + (y i /2 + / i i )(n i ), 



where are restricted by the stability requirement of 
the order parameters. All expectation values in @ and 
may be straightforwardly calculated according to the 
variational wave function IjlQI) . 

The phase diagram may obtained by solving the con- 
straint equations for the stable order parameter and min- 
imizing the variational energy To save the space, 
we do not discuss strongly asymmetric case with t-\ /V\ <C 
t[/Vi in which one-component bosons may in the Mott 
phase while the other is in the superfluid phase. We fo- 
cus on the species-symmetric model with t-\ = t± = t 
and Vf = Vj. = V. The symmetry requires that the 
order parameters are species-independent and leads to 
/.tj = hi = \x and \ — C = 7r/2. The variational param- 
eters are reduced to 9 and £ only. Corresponding to the 
stability of non-zero (a-f) = (aj.), (a?) = (of), or (a^ai), 
the constraint identities are 



{zt + n+ U/2) ( 1 + sin - + cos - J = U (sin - + cos - 
2 2 J V 2 2 



zt sin 2 - 1 + sin - + cos - + 2(ji + U/2) cos - = 



sin 2 —11 + sin — + cos - J + 2fi sin - = 0. (12) 
— \ — — / — 



E 



tz 



1 9 0/ 9 £ 9 £ 

- cos 2 - V cos 2 - + U sin 2 - 

2 2\ 2 2 



(13) 
V 



Two cases are analytically solvable. One is the hard-core 
limit in which V/U — * oo and another is the critical point 
of the phase separation, U — V. In the hard-core limit, 
the energies are minimized at £ = tt. This gives (a 2 ) = 
as it should be in the hard-core limit. The atomic super- 
fluid phase (AS/MS) corresponding to the first constraint 
in i|12[l is confined — 1 < U/Atz < 1 while the molecular 
superfluid phase (MS/ AS) to the last constraint in 112|) 
is confined < U/4tz < 2. The variational energies cor- 
responding to these two order parameters in the whole 
regime of the interaction U may be given by in a similar 
form as (JJJ . Comparing these two energies for any given 
U/4zt, the phase diagram in the hard-core limit for the 
symmetric model can be determined (see Fig. 2(a)). The 
four phases can detected: the molecular superfluid phase 
and the mixture of the stable molecular superfluid and 
unstable atomic superfluid phases, these two phases cor- 
responds to the s-wave superconductivity. The mixture 
of the stable atomic superfluid and unstable molecular 
superfluid phases, which corresponds to the preparing 
(but (i-wave) for the electrons in high T c materials ; and 
the Mott insulator phase. The boundary of the second 
and third ones is in U/Azt = 0.5. 

For any finite U/V , to keep the variational energy hav- 
ing a lower bound, the second and third constraints in 
(|12fl can only be satisfied if 9 = 0, and /i = —U/2 and 
0, respectively. In the critical point of phase separation, 
U = V , E2 = iV/2 for the second constraint and E3 = V 
for the third one and then the dimerized phase is always 
unfavorable. The energy corresponding to the first con- 
straint is minimized at £ = n/2. We compare £3 and E\ 
which are in a similar form as Q and the phase diagram 
can be depicted (see Fig. 2(b)). 

For a general finite U/V < 1, the variational energies 
may be numerically minimized and the phase diagrams 
are plotted in Fig. 2(b) for U > 0. It is seen that the 
'MS/AS-D' phase can only exist in the hard-core limit. 
For a finite V/U, this phase is suppressed and turns to 
the 'MS' phase. However, as V reduces, the atomic su- 
perfluid phase 'AS/MS-D' grows in both directions. In 
our variational state, the suppression of the Mott phase 
is not so strong when V/U < 1.73. 

We now discuss the asymmetric case in the hard-core 
limit. For t-\ 7^ t±, we find that the atomic superfluid 



state with a stable (an) exists for 



At tan xo 



< 



U 



< 



2 sin xo 
1 



2z(t T +tj.) ^ 2 

of the equation 4(1+sinx) 

— - At tan x and At 



2sinxo 2 

Att ^ nx ° , with xo is the solution 

U sin x 1 A£ cos x 



2z(t r +t l )(l+sin X ) 



2(l+sin X ) ^ 

Solving the equa- 



z_(*T+*l) —-a— *T+*i' 

tion of x, one finds that there is only a minor modifica- 
tion to that uses the simple variational method taking 
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H„ = 0. The latter gives \U\/2z{t^ + t t ) < 1 9] and the 
former, even for At = ±1, \U\/2z(t ] +t{) < 0.97. At the 
other critical points, one also finds the critical values of 
\U\/2z(t^ + ti) does not leave far from their symmetric 
ones. Thus, we can think the basic feature of the phase 
structure of superfluid phase in the hard-core limit can 
mostly be described by the symmetric model. For finite 
V a , the situation may be more complicated, e.g., for the 
strongly asymmetric case, there is a mixture between the 
Mott insulator in the one species and the superfluid in 
the other (see, 0, Chen and Wu at 0). We shall leave 
this for coming works. 

Finally, we comment on the counterflow superfluid (see 



Kuklov et al in Q ) • The corresponding order parameter 
is (ajaj.) = 5 [sin 2 | sinx + a/2 cos 2 | sin £ cos 2 f (sin | + 

cos|)]. For 6 = 7r, the Mott state, (atdj.) = (1/2) sin %. 
However, this mean field order parameter is not stable 

because the Heisenberg equation of motion — i((ajajj) = 

V(((rif — njjaja^)) = V^(ajaj) for the Mott state. In 
the symmetric model, it holds for arbitrary U and V0. 
Thus, from a mean field point of view, the counterflow 
superfluid is stable only in a time scale t <C 1/V. 

This work was supported in part by the National Sci- 
ence Foundation of China. 
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Fig. 1 (a) The variational energies associated with differ- 
ent stable order parameters for the one-component Bose- 
Hubbard model, (b) The corresponding phase diagram. 
'D' is the dimerized phase with (a 2 ) = 1/2 and (a) = 0; 
'SF-D' the superfluid phase with a stable (a) 7^ and un- 
stable (a 2 ) ^ 0; Mott' is the Mott insulator phase with 
(a) = {a*)=0. 

Fig. 2 The phase diagrams for different V/U. MS' 
denotes the molecular superfluid phase with a stable 
(aj&x) = i and other order parameters vanish. 'MS/AS- 
D' denotes all order parameters non- vanishing but only 
(ffll-fflj.) stable. 'AS/MS-D' does similarly but (a CT ) stable, 
(a) The hard-core limit, (b) Finite V/U. 
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